Direct Measurement of Supernova Neutrino Emission Parameters 
with a Gadolinium-Enhanced Super-Kamiokande Detector 
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The time-integrated luminosity and average energy ol the neutrino emission spectrum are essential 
diagnostics of core-collapse supernovae. The SN 1987A electron antineutrino observations by the 
Kamiokande-II and 1MB detectors are only roughly consistent with each other and theory. Using 
new measurements of the star formation rate history, we reinterpret the Super-Kamiokande upper 
bound on the electron antineutrino flux from all past supernovae as an excluded region in neutrino 
emission parameter space. A gadolinium-enhanced Super-Kamiokande should be able to jointly 
measure these parameters, and a future megaton-scale detector would enable precision studies. 
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When a massive star dies, its core collapses and re- 
bounds, producing an outgoing shock wave that should 
eject the stellar envelope, causing the optical supernova, 
and leaving behind a neutron star remnant. However, in 
simulations, the shock wave stalls, leading to the whole 
star collapsing into a black hole, failing to produce an 
optical supernova or spread its heavy-element yields 0. 
Since the required explosion energy is only ~ 1% of the 
emergent neutrino energy, a full accounting of the neu- 
trino emission is essential for understanding supernovae. 
Further, in the Bethe- Wilson delayed explosion model, 
the neutrinos revive the shock 0. Resolution of the su- 
pernova problem would also have profound implications 
for the history of stellar evolution and nucleosynthesis. 

The weak interactions of neutrinos, which allow them 
to reveal the dynamics deep within the exploding star, 
also make their detection challenging. The last nearby 
supernova, SN 1987A, occurred in the Large Magellanic 
Cloud at 50 kpc, and ~ 20 neutrinos were detected Q 
preceding the optical supernova, confirming our basic un- 
derstanding of the explosion 4]. However, even taking 
into account the small statistics, the fitted ranges for the 
time-integrated luminosity and average energy are per- 
plexing, showing clear discrepancies among the experi- 
mental detections and theory [H [|| . A Milky Way super- 
nova would yield many events in present detectors, but 
the expected supernova rate is only ~ 3 per century. We 
have shown that with proposed megaton-scale detectors, 
it will be possible to build up the spectrum by detecting 
neutrinos one or two at a time from supernovae within 
10 Mpc, at a rate as large as ~ 1 neutrino per year 0. 

Here we propose a new approach, which could begin 
immediately, if the existing Super-Kamiokande (SK) de- 
tector were modified by the addition of gadolinium to 
greatly reduce backgrounds, as proposed by Beacom and 
Vagins H, |j| . We consider the spectrum of the Diffuse 
Supernova Neutrino Background (DSNB) H3, EH E3, 
as the observable. The DSNB predictions depend on the 
redshift evolution of the supernova rate, which is sepa- 
rately measurable and increasingly well known, and the 
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FIG. 1: (Color online) Joint limits on the time-integrated 
luminosity L„ and the spectrum average energy (E„) for elec- 
tron antineutrinos. The two contours are the allowed regions 
at 90% C.L. from the SN 1987A analysis of Ref. Q and our 
shaded region corresponds to the SK 90% C.L. upper limit on 
the DSNB flux 16]. 



neutrino emission per supernova, the object of our study. 
While the received neutrino spectrum will be redshiftcd, 
it will have relatively high statistics, up to several events 
per year in SK. Recently, the DSNB uncertainties from 
the star formation rate history |l4l Il5| narrowed enough 
that it is now sensible to reinterpret the SK flux limit [lJ| 
as an exclusion region in the plane of the time-integrated 
luminosity and average energy, which can be directly 
compared to the allowed regions from SN 1987A. Antic- 
ipating further improvements in the astronomical data, 
we show that a gadolinium-enhanced SK should be able 
to usefully constrain the emission parameters in much of 
the interesting range. 
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Supernova 1987A Signal. — One of the triumphs of 
astrophysics, nuclear physics and particle physics was 
the detection of neutrinos from SN 1987A, so far the 
only astrophysical source besides the Sun seen with neu- 
trinos. The ~ 20 events in the Kamiokande-II (Kam- 
II) and 1MB detectors are assumed to be mostly in- 
verse beta events, v e + p — > e + + n; the cross sec- 
tion a ~ -Ey, with the positron carrying nearly the full 
neutrino energy 01 ■ The number of detected events 
N det ~ {LJ{E U )) ■ (E v ) 2 ~ L v ■ (E u ), where L v is the 
time-integrated luminosity of the electron antineutrinos 
and {E„) is the average energy of the neutrino emission 
spectrum. The average detected energy (Edet) ~ {Eu)- 
The combination of these two constraints explains the 
banana-shaped allowed regions shown in Fig. ^ taken 
from the full spectrum analysis of Ref. [5j (those authors 
assumed a Maxwell-Boltzmann thermal emission spec- 
trum). We show only the 90% C.L., which is appropriate 
for this level of precision, to avoid cluttering the figures. 

At least three puzzling features of the SN 1987A data 
still stand out. First, the fits to the Kam-II and 1MB data 
for the neutrino emission parameters barely overlap, due 
to the disagreement on the spectra 0, |(| . Second, the 
results disagree with the canonical expectations [a, @j 
conservatively indicated by point D in Fig. ^ This cor- 
responds to a neutron star binding energy of 3 x 10 53 erg, 
assumed shared equally among the six flavors, and an ef- 
fective received v e temperature of about 5 MeV. Both 
the Kam-II and 1MB data allow very high luminosities, 
perhaps reflecting a larger neutron star binding energy 
and/or a violation of its assumed equipartition among 
flavors. Both, but especially Kam-II, allow very low av- 
erage energies, especially if neutrino mixing with higher- 
temperature flavors is taken into account (i.e., with tem- 
peratures possibly as large as 8 MeV). Third, both the 
Kam-II and 1MB results are in significant disagreement 
with model- independent tests of the angular distributions 
of the detected events [TllTl l2Cj . 

We emphasize that the Kam-II and 1MB results are 
roughly consistent with each other and theory; still, there 
are puzzling issues raised which cannot be answered with- 
out new data. Also, all supernovae may not be alike, and 
the DSNB results will reveal the average neutrino emis- 
sion parameters of supernovae, possibly more relevant for 
cosmological applications. We therefore stress that the 
adoption of SN 1987A as a template for DSNB studies is 
undesirable. 

SK DSNB Limit. — To predict the DSNB flux, one 
needs only the core-collapse supernova rate as a func- 
tion of redshift z, convolved with the neutrino emis- 
sion per su per nova, taking into account redshift ef- 
fects [T3, Ull [12; 03 ■ At present, the former is calculated 
from the measured star formation rate and the stellar ini- 
tial mass function, which determines the fraction of stars 
that end their lives as core-collapse supernovae. We base 
our results primarily on the GALEX star formation rate 
data [T^ | , for which the normalization uncertainty is now 
at the ~ 30% level. This yields results similar to those 
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FIG. 2: (Color online) DSNB detection spectra for selected 
parameters (solid and dashed curves), efficiency-corrected SK 
data (points with error bars), and detector background (solid 
steps), all in counts per 4 MeV, per year, per SK fiducial 
volume. The values for the representative points A, B, C 
and D are given in Table Q] SK is so far sensitive only to the 
dark shaded region above 18 MeV due to high backgrounds at 
lower energies (not shown). With the addition of gadolinium, 
these backgrounds in the range 10-18 MeV would be removed, 
and that shown reduced by a factor ~ 5, opening up also the 
light shaded region for analysis 0,0. 



of the Concordance Model of Ref. 0] , which was shown 
to be consistent with the latest measured star formation, 
thermonuclear (type la) supernova, and core-collapse su- 
pernova (types II, lb, and Ic) rates, as well as other data, 
and which predicts a DSNB flux just below the present 
SK limit. (See also the more recent Ref. [3.) The DSNB 
prediction thus depends first on a purely astronomical 
factor, which is already measured well and will ultimately 
be precisely and unambiguously measured through direct 
data on the supernova rate versus redshift (note that op- 
tically failed supernovae with substantial neutrino emis- 
sion would increase the core-collapse rate) [l3| . The sec- 
ond factor, the average neutrino emission per supernova, 
must be measured directly, and this is our focus. 

In 2002, SK reported their DSNB flux upper limit for 
electron antineutrinos with a detection energy thresh- 
old of 18 MeV as 1.2 cm" 2 s" 1 (90% C.L) through 
non-detection of excess counts above background fluctua- 
tions ^(| . Due to the rising background and falling signal 
with increasing energy, almost all of the statistical power 
derives from the first two bins in Fig. [21 which is why 
the SK flux limit was the same for DSNB models with 
different spectral shapes. With the present statistics, it 
is enough to use just these two bins to limit the signal, 
noting that the other bins fix the background normaliza- 
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TABLE I: The values for the points A, B, C (near the SK 
upper bound) and D (canonical values) of the figures. 



Point 
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L v [10 52 erg] 


Sensitivity 
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15 
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25 
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15 
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tion. While at the time of the SK analysis, the DSNB 
models differed significantly in their normalization, the 
latest astronomical data greatly restricts this freedom, 
and will eventually eliminate it, modulo the differences 
in neutrino emission per supernova that we want to test. 
Thus it now makes sense to reinterpret the SK event rate 
limit (~ 2 yr _1 in 18-26 MeV in terms of the su- 
pernova electron antineutrino emission parameters, the 
time-integrated luminosity and average energy. The re- 
sult of our analysis is that the shaded region in Fig. 
is excluded. (We assumed a thermal emission spectrum 
of the form used in Ref. |6j, taking a — 3, which corre- 
sponds to a somewhat "pinched" spectrum.) This does 
not yet reach the allowed regions deduced from the SN 
1987A data, but it is encouragingly close. Neutrino mix- 
ing can blend the initial D e spectrum with higher-energy 
v^jvr spectra. We are limiting an effective composite 
spectrum, which will be dominated by the harder spec- 
trum Thus our analysis is conservative, in that with 
neutrino mixing, the interpretation of the DSNB bound 
would be more constraining (e.g., Ref. [l3j]). 

Beyond the two supernova neutrino emission parame- 
ters used here, there is also the question of the spectrum 
shape, and whether it is distorted from thermal by being 
"pinched" or "anti-pinched" jaL The SK energy thresh- 
old of 18 MeV is used in Ref. is high, especially not- 
ing that redshifts z < 1 are relevant; the SK limit is thus 
based on energies ~ 20 — 40 MeV in emission, where the 
spectrum uncertainties are largest. Indeed, this is part 
of the motivation for lowering the SK energy threshold, 
so that the detected events would correspond to emission 
from the better-understood spectrum peak region. 

One can get more insight by examining three points, 
A, B, and C, shown in Fig. ^ and Table [I] which are 
at the edge of detectability. We also consider a point 
D, which is often regarded as the canonical values for j/ e 
emission before neutrino mixing. The DSNB spectra for 
these points are shown in Fig. |3 together with the SK 
data and background expectations (dominantly from the 
decays of sub-Cerenkov muons produced by atmospheric 
neutrinos). The three points A, B, and C correspond 
to almost equal yields (comparable to the fluctuations 
in the backgrounds) in 18-26 MeV, where the present 
SK sensitivity is greatest. The point D produces fewer 
signal events and is safely allowed. Note that these spec- 
tra are quite different at lower energies. Thus it is clear 
that in order to make the necessary progress over the 
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FIG. 3: (Color online) Possible 90% C.L. measurements of 
the emission parameters of supernova electron antineutrino 
emission after 5 years running of a gadolinium-enhanced SK 
detector. 



present background-limited search, SK must reduce both 
the background rates and the energy threshold. 

The points A, B, and C, besides representing three 
very different possibilities, are also the most favorable for 
being probed by SK in the near term. (Below, we discuss 
the outlook for models with a lower neutrino emission per 
supernova.) These points are chosen just with respect to 
what is allowed by the DSNB, deliberately not taking 
into account other possible constraints, so that the im- 
pact of our results can be clearly seen. Additionally, it is 
important to keep an open mind about what the true pa- 
rameters are, given that (a) numerical supernova models 
fail to explode pj, and (b) it remains possible that SN 
1987A was very different from an average supernova. 

Points A and C are relatively extreme, compared to 
the more canonical point D. However, possibilities like A 
might occur if the neutron star binding energy and the 
fraction of this energy carried away by v e (i.e., a violation 
of the usually assumed equipartition of the energy among 
the six flavors) are both larger than expected [l|; more 
directly, A is very close to part of the 1MB allowed region 
for the SN 1987A data. Possibilities like C correspond to 
a temperature of about 8 MeV, which is within the range 
considered in many papers, especially if neutrino mixing 
is taken into account. 

Gadolinium-Enhanced SK Sensitivity. — In SK, 
DSNB z7 e would be detected by the inverse beta reac- 
tion v e + p — > e+ + n on free protons 0. At present, 
this is a (positron) singles search, for which there are very 
large background rates ^(J- With dissolved gadolinium, 
SK could also detect the neutron via its radiative cap- 
ture (TiEj ~ 8 MeV) , which would give a tight temporal 
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and spatial coincidence for the signal events |8|, |9j . This 
would reduce the sub-Cerenkov muon decay background 
shown in Fig. 0by a factor ~ 5, and would remove the 
spallation backgrounds in the range 10-18 MeV; below 
about 10 MeV, the reactor v e signal suddenly becomes 
dominant 0, Ej . We emphasize that SK would be work- 
ing well within its design range when detecting these 
positrons and neutron captures, i.e., at these energies 
the detection efficiency is very high, nearly constant, and 
well-measured through calibrations. In contrast, in the 
detection of neutrinos from SN 1987A, both Kam-II and 
especially 1MB had events detected at energies where the 
efficiency was low and/or varying quickly |3|]. 

Besides increasing the signal rate and improving the 
ability to test the DSNB spectral shape, the capability 
for neutron tagging could allow a rate-limited, instead of 
background-limited, search, so that the sensitivity could 
improve linearly with detector exposure. To examine the 
prospects for a gadolinium-enhanced SK, we again con- 
sider the four points in Table [I] For each, we simulate 
the expected neutrino spectra over a 5-year period (i.e., 
5 times the yields shown in Fig. [21 noting the background 
reduction). We fit the spectra of DSNB and background 
events simultaneously, and compute the x 2 in 4-MeV 
bins, as in Ref. 16]. With these statistics, the Gaussian 
X 2 is adequate to draw the 90% C.L. contours. 

In Fig. |3| we show the expected determinations of 
the physical parameters at 90% C.L. The spectra cor- 
responding to A, B, and C, while presently indistin- 
guishable above 18 MeV, would be clearly separable in 
a gadolinium-enhanced SK (see Fig. |2J). With a lower 
energy threshold, there is greater sensitivity to both the 
time-integrated luminosity and average energy through 
the spectral shape as well as normalization. At differ- 
ent points, the relative sensitivity to the two parameters 
changes, as listed in Table [i] and shown in Figs. and [3] 
If we allowed the supernova rate history to be free, then 
the normalization of the supernova rate would be de- 
generate with the time-integrated luminosity, while the 
rate of increase with redshift would be degenerate with 
the average energy. Thus near point A, inaccuracies in 
the normalization would have little impact, and likewise 
near point C for inaccuracies in the rate of increase with 
redshift. While we think that the supernova rate uncer- 
tainties will play a minor role compared to those on the 
supernova emission parameters, this may make the de- 
termination of one variable more robust than the other. 

If the true parameters, even taking neutrino mixing 
into account, are closer to the canonical values (point 
D), then the detection rate will be lower, and hence the 
allowed region larger. (From the size and shape of this 
contour relative to the others, one can estimate how the 
contours for other points would look.) Even in the case 
of point D, in 5 years a gadolinium-enhanced SK would 
detect ~ 15 events, comparable to the SN 1987A yield, 
which should help distinguish between the Kam-II and 
1MB solutions. If the Kam-II region is correct, then the 
detection rate will be even lower; note that an outcome 



that disfavored the 1MB region might be viewed as se- 
lecting the Kam-II region, given the prior information 
from SN 1987A. If the supernova neutrino emission pa- 
rameters are both low, a megaton-scale detector may be 
necessary to accumulate good statistics. 

Discussion. — The ~ 20 neutrinos from SN 1987A 
provided a rough confirmation of core-collapse supernova 
neutrino emission, and hence of the dynamics of the ex- 
ploding star in the first several seconds after collapse. 
These details of core-collapse events will forever remain 
invisible with photons, but can be revealed by neutrinos, 
if they can be detected. While very challenging, it is hard 
to overstate the importance of this goal. Since nearly 20 
years after SN 1987A, we have no further direct informa- 
tion on supernova neutrino emission, techniques besides 
waiting for a Milky Way supernova must be considered. 

We propose that with the rapidly improving preci- 
sion of the astronomical data, it will be possible to use 
measurements of the DSNB to constrain the v e emission 
per supernova. This information will be limited, in that 
neither the dynamic timescales nor the emission in the 
other neutrino flavors can be measured; only a Milky 
Way supernova can provide those data. However, even 
with just v e , the time-integrated luminosity and aver- 
age energy will constrain the explosion energy and proto- 
neutron star opacity, especially if reasonable assumptions 
are made about the other flavors. Despite these caveats, 
and the limited statistics, we stress that this technique is 
unique in that in a very short time it could begin provid- 
ing steadily better clues to the mysteries of SN 1987A. 

Our results in Fig.[3]show that a gadolinium-enhanced 
SK detector would have useful sensitivity to an interest- 
ing range of supernova emission parameters. The recog- 
nition that the astrophysical uncertainties are small and 
quickly diminishing allowed us to reinterpret the SK flux 
constraint [16j in terms of the neutrino emission per su- 
pernova. As a practical matter, we encourage SK to rep- 
resent future results in this way, as it is more directly 
connected to the measured event spectrum than the inte- 
gral of the flux above a given energy (in particular, since 
the latter does not contain the important weighting by 
the detection cross section), i.e., is less model-dependent. 

If a Milky Way supernova is detected, this would in- 
crease the value of the proposed DSNB measurement. 
The comparison of results could probe whether the neu- 
trino emission from core-collapse supernovae is as uni- 
form as presently assumed. Alternatively, it could test 
the measured core-collapse rate history [l2j , and whether 
there is an additional neutrino background from explo- 
sions which fail |l3j . emitting neutrinos but not creating 
an optical supernova, just as is seen in simulations Q. 

Proposed megaton-scale detectors would greatly ex- 
tend the sensitivity to these and more general spectra, 
and could bring precision to the measurement, due to 
the very high DSNB statistics. Such detectors could also 
allow the accumulation of events from identified super- 
novae within 10 Mpc 0. While the statistics of the latter 
would be relatively lower, that spectrum would not be 
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redshifted, nor dependent on the evolution of the cosmic 
supernova rate, and hence would be complementary to 
the DSNB spectrum. 

The neutrino emission per supernova is also impor- 
tant for understanding nucleosynthesis, especially of the 
heavy elements beyond iron, which are believed to be 
formed only in core-collapse supernovae, and which re- 
quire special conditions that may be importantly affected 
by the neutrinos |2l|. In addition, Yoshida et al. [z| have 
recently shown that the yield of the light element n B 
constrains the neutrino emission parameters to be close 
to the canonical values, which is favorable for confirma- 
tion by direct detection. Combining the nucleosynthesis 
results [12 with future sensitivity to the DSNB electron 
antineutrino flux (as stressed here) E | fToL fTTL H3 . . the 
DSNB electron neutrino flux [lajl9|, and the summed 
spectrum of nearby supernovae [7J will provide comple- 
mentary and restrictive probes of the details of supernova 
neutrino emission and the history of stellar birth, life, and 
death. 



Acknowledgments. — We thank L. Strigari and 
M. Vagins for discussions. This work was supported by 
The Ohio State University and NSF CAREER grant No. 
PHY-0547102 to J.F.B.; S.A. was also supported by the 
Japan Society for the Promotion of Science. 



[1] M. T. Keil, G. G. Raffelt and H. T. J anka, Astrophys. 
J. 590, 971 (2003) astro-ph/0208035 ; T. A. Thomp- 
son, A. Burrows and P. A. Pinto, Astrophys. J. 592, 434 
(2003) |astro-ph/0211194| ; M. Liebendoerfer, M. Rampp, 
H. T. Janka and A. Mezzacappa, Astrophys. J. 620, 840 
(2005) astro-ph/0310662 . 

[2] H. A. Bethe and J. R. Wilson, Astrophys. J. 295, 14 

(1985) . 

[3] K. Hirata et al, Phys. Rev. Lett. 58, 1490 (1987); 

R. M. Bionta et al, Phys. Rev. Lett. 58, 1494 (1987). 
[4] W. D. Arnett, J. N. Bahcall, R. P. Kirshner and [12] 

S. E. Woosley, Ann. Rev. Astron. Astrophys. 27, 629 

(1989); D. N. Schramm and J. W. Tmran, Phys. Rept. [13] 

189, 89 (1990). 

[5] B. Jegerlehner, F. Neubig and G. Raffelt, Phys. Rev. D [14] 
54, 1194 (1996) astro-ph/9601111. 

[6] A. Mir izzi and G. G. Raff elt, Phys. Rev. D 72, 063001 [15] 
(2005) astro-ph/0508612 . [16] 

[7] S. Ando, J. F. Beacom and H. Yiiksel, Phys. Rev. Lett. 

95, 171101 (2005) astro-ph/0503321 . [17] 

[8] J. F. Beacom and M. R. Vagins, Phys. Rev. Lett. 93, 

171101 (2004). hep-ph/0309300 . [18] 

[9] M. R. Vagins (Principal Investigator) , Department of En- 
ergy 2003 and 2005 Advanced Detector Research pro- [19] 
gram grants; M. R. Vagins, talk at "NNN05: Next Gen- 
eration of Nucleon Decay and Neutrino Detectors," 7-9 
April 2005, Aussois, France, http://nnn05.in2p3.fr/ [20] 
[10] S. And o, K. Sato and T. To tani, Astropart. Phys. 18, 307 

(2003) astro-ph/0202450 ; S. Ando and K. Sato, Phys. [21] 
Lett. B 559, 113 (2003) astro-ph/0210502 . 
[11] G. S. Bisnovatyi-Kogan and Z. F. Seidov, Sov. As- 
tron. 26, 132 (1982); L. M. Krauss, S. L. Glashow and 
D. N. Schramm, Nature 310, 191 (1984); S. E. Woosley, 
J. R. Wilson, and R. Mayle, Astrophys. J. 302, 19 [22] 

(1986) ; R. A. Malaney, Astropart. Phys. 7, 125 (1997) 
|astro-ph/9612012| ; D. H. Hartmann and S. E. Woosley, 
Astropart. Phys. 7. 137 (1997); M. Kaplinghat, 



G. Steigman a nd T. P. Walker, Phys. Rev. D 62, 
043001 (2000) astro-ph/9912391 ; M. Fukugita and 
M. Ka wasaki, Mon. Not . Roy. Astron. Soc. 340, L7 

(2003) astro-ph/0204376 ; L. E. Strigari, M. Kapling- 
hat, G. Steigman and T. P. Walker, JCAP 0403, 007 

(2004) astro-ph/0312346 ; S. Ando and K. Sato, New 
J. Phys. 6, 170 (2004) astro-ph/0410061 ; C. Lunardini, 
|astro-p h/0509233l F. Daigne, K. A. Olive, P. Sandick 
and E. Vangioni Phys. Rev. D 72, 103007 (2005) 
astro-ph/0509404 . 

S. Ando, Ast rophys. J. 607, 20 (2004) 
astro-ph/0401531 . 

L. E. Strigari, J. F. Beac om, T. P. W alker and P. Zhang, 
JCAP 0504, 017 (2005) [astrb -ph/0502150 . 
D. Schiminovich et al. Astrophys. .1. 619. L47 (2005) 
astro-ph/0411424 . 



A. M. Hopkins and J. F. Beacom, astro-ph/0601463 
M. Malek et al, Phys. Rev. Lett. 90, 061101 (2003) 
hep-ex/0209028 . 

P. Vog el ancHL F. Bea com, Phys. Rev. D 60, 053003 
(1999) hep-ph/9903554 . 

J. F. B eacom and L. E. S trigari, Phys. Rev. C 73, 035807 
(2006) hep-ph/0508202 . 

A. G. Cocco, A. Ereditato, G. Fiorillo, G. Mangano and 
V. Pettorino, JCAP 0412, 002 (2004) hep-ph/0408031 ; 
C. Lunardini, hep-ph/0601054 

M. L. Costantini, A. Ianni and F. Vissani, Phys. Rev. D 
70, 043006 (2004) [astro-ph/0403436 . 
J. Fetter, G. C. McLaughlin, A. B. Balantekin 
and G. M. Full er, Astropart. Phys. 18, 433 (2003) 
hep-ph/0205029 ]; Y. Z Qian, Pr og. Part. Nucl. Phys. 
50, 153 (2003) |astro- ph/0301422]; A B. Balantekin and 
H. Yiiksel, New J. Phys. 7, 51 (2005) astro-ph/0411159 . 
T. Yoshida, T. Kajino and D. H. Hartmann, Phys. Rev. 
Lett. 94, 231101 (2005) astro-ph/0505043 . 



